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Abstract 

An analytical nonadiabatic approach has been developed to study the dimerization 
gap and the optical absorption coefficient of the Su-Schrieffer-Heeger model where the 
electrons interact with dispersive quantum phonons. By investigating quantitatively the 
effects of quantum phonon fluctuations on the gap order and the optical responses in this 
system, we show that the dimerization gap is much more reduced by the quantum lattice 
fluctuations than the optical absorption coefficient is. The calculated optical absorption 
coefficient and the density of states do not have the inverse-square-root singularity, but 
have a peak above the gap edge and there exist a significant tail below the peak. The 
peak of optical absorption spectrum is not directly corresponding to the dimerized gap. 
Our results of the optical absorption coefficient agree well with those of the experiments 
in both the shape and the peak position of the optical absorption spectrum. 
PACS number(s): 7L38.-k; 7L20.Rv; 78.20.Ci; 7L45.Lr 
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I Introduction 



Since the pioneering work of Su, Schrieffer and Heeger (SSH) the SSH model has been 
widely used to investigate the properties of conducting polymers, especially polyacetylene 
[(CH) X ] [|J. At half-filling, this system, as many quasi-one-dimensional materials, exhibits 
an instability against a periodic lattice distortion due to the Peierls dynamics. Because 
of a twofold degenerate ground state, this system allows for the soliton excitation with an 
associated electronic state at the center of the energy gap [||, While considerable effort has 
been made to understand excitation energies of solitons and polarons in these compounds, it 
is theoretically and experimentally significant to develop an accepted theory for their optical 
absorption spectrum. The recent exploitations of this model in studying the properties of 
nanotubes [||, ^, f?| renew further interesting of theoretical studies on the intrinsic physics 
of this model. In the previous works on conducting polymers, along with the expanding 
of measurements, theoretical attempts have been made to study the energy gap Q, low- 
lying excitation ||, [|, ||, and optical response || |l(J of this system. However, there exist 
discrepancies between theories and experiments, (i) Within the adiabatic approximation, 
the calculated optical absorption coefficient of the perfectly dimerized lattice has an inverse- 
square-root edge singularity at u = 2A and there is no absorption inside the gap. But the 
observed optical absorption spectrum of polyacetylene is quite different jll], 12, The 
singularity is absent, and there is a significant tail below the maximum, (ii) The energy gap 
deduced from the absorption edge is smaller than the activation energy of the dc conductivity 
0,0- 

In (CH) X the zero-point motion of the lattice Su (=0.03 is comparable to the ampli- 
tude of the lattice Peierls distortion uq (=0.03 or 0.035 ||l7| ). Although it has been stated 
that the quantum lattice fluctuations cause a slightly small reduction of the dimerization or- 
der parameter JO], [l^, |2(J, if 5u ~ uq, as was pointed out in Ref. [Q, ^5| , one might wonder 
why there should be a clearly defined gap in the spectrum at all and expect the optical absorp- 
tion spectrum to peak at the correspondingly renormalized gap edge. In nonadiabatic case, 
the phonon frequency is finite. Generally speaking, the nonadiabatic effect will suppress 
the state order parameters of the system (T^j. As the density-of-states (DOS) of electrons is 
concerned, the results of adiabatic approximations also have inverse-square-root singularity 
at the gap edge. By considering the nonadiabatic effect, the singularity may disappear Q]. 
The influence of the phonon frequency on DOS, the optical-absorption spectrum, and other 
order parameters in the range from u> n = to uj n — > oo should be studied for understanding 
the physics of electron-phonon interactions in nonadiabatic case. 

The theoretical analysis of the coupled electron-phonon system for all values of the ionic 
mass and the electron-phonon coupling constant is a very difficult problem, especially when 
the quantum lattice fluctuation and the lattice dynamic (quantum phonon) effects are taken 
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into account. This problem has been investigated by using the Monte Carlo simulation 
Jl8|, ni!J, perturbation calculation p3], Green's- function technique [fRI], renormalization group 



analysis p3| , |24}| , variational method of the squeezed-polaron wave-function [25], etc. In 
a resent work [26], this problem was studied by considering the spinless Holstein model 
and obtained some interesting results. In this paper, by considering SSH model, the more 
realistic case with spin-i electrons, we focus on the properties of the effects of quantum 
lattice fluctuations on the dimerization energy gap and the optical-responses of the system 
with the view of understanding the measurement evidences and of interpreting why and how 
the dimerization energy gap is smaller than the optical absorption gap. Our analyses show 
that at finite phonon frequency, the peak position of optical absorption spectrum is not 
directly corresponding to the dimerized gap. The dimerization gap is much more reduced 
by the quantum lattice fluctuations than the optical absorption coefficient is. The effects of 
quantum lattice fluctuations on the dimerization gap and on the optical gap are essentially 
different, especially when phonon frequency uj n is small. The reason for this is discussed 
and the relationship between the peak position of optical absorption spectrum and the true 
dimerization gap is obtained quantitatively. 



II Model 

We start from the SSH model ill 



H = ~ 53 + a ( u l ~ u l+l)K c l s c l+l,s + 4+1 s C hs) + 53 



Ls I 



(1) 



The notations in it are as usual |2j (throughout this paper, we set fi = ks = 1). 

In Hamiltonian ([[]) the operators of lattice modes, ui and pi, can be expanded by the 
phonon creation and annihilation operators b\ and b q , and after Fourier transformation to 
momentum space, the H becomes 

h = 53^ c L c m + 53^ + 2) + -f= 53 0(fc+?>*O(&i« + m4 +<? , s c m> ( 2 ) 

where = —2tcosk is the bare band structure, TV" is the total number of sites. The phonon 
frequency 

u q = y(4i^/M)sin 2 |, (3) 

and the coupling function 

g(k + q,k) = i2a^l/(2Muj q )[sm(k + q) - sink}. (4) 

In the SSH model, the electron couples to the difference between the phonon amplitudes on 
the two neighboring sites and the form of the lattice vibration energy leads the phonons to 
be dispersive. 
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Among the models for one-dimensional systems the Holstein model [27] and SSH model 
are the two typical electron-phonon coupling Hamiltonian studied by many previous authors. 
The Holstein model is for the on-site electron-phonon interaction with dispersionless phonons 
coupled with electron density operator and has played a central role in the polaron problem, 
while the SSH model is for the on-bond electron-phonon interaction and has been successful 
in describing the essential physics of conducting polymers. Phonon dispersion and different 
structures of the coupling terms lead to important differences in the physics of the two 



models. For the Holstein model, although numerical simulation [28| has been performed in 
confirmation of the supposition that it is likely the Peierls gap is more reduced by the quantum 
lattice fluctuations than the dimerization or the charge-density-wave amplitude, in other 
words, the quantum lattice fluctuations have a much strong effect on the Peierls gap than 
on the amplitude of the Peierls distortion, due to the difficulty in preparing practicably pure 
sample of molecular crystal materials, few of observation data are available for the comparison 
with theoretical results. Other models somewhat intermediate between the SSH and Holstein 



models, such as the model for MX chain materials [29] and the model for electron- libron 
coupling in Polyaniline J3C|], because of their complex forms with Coulomb interactions and 
steric potential accounting for the ring torsion angle, will make the presentation of our analysis 
to be more prolix and inexplicit. The enormous measurement data of optical absorption 
spectrum in conducting polymers and the relatively simple form of SSH model render it the 
favorable system for our investigation into the nonadiabatic effect on the dimerization gap 
and the optical absorption spectrum of the electron-phonon coupling systems. The previous 
researches on this problem in SSH model by using a static Gaussian-random potential method 
to model the lattice disorder field [15, 31] and a configurational-coordinate method to treat the 
anharmonic lattice fluctuations [32] were phenomenological theories and have not obtained 
the true Peierls dimerization gap. It is the purpose of this paper to develop an nonadiabatic 
approach to study analytically this problem simply from the model Hamiltonian. 



Ill Effective Hamiltonian 

The electron-phonon coupling and the Peierls dimerization are two main respects in our 
studying on the SSH model. First, in order to take into account the electron-phonon coupling, 
an unitary transformation is applied to H, 

H' = exp(S)Hexp(-S), (5) 

where the generator S is 

s = 1 ^ 9(k + q,k) 6{k + ft k){ ^ _ bq) j (6) 

q,k,s U 1 
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Here, we introduce a function S(k',k) which is a function of the energies of incoming and 
outgoing electrons in the electron-phonon scattering process and it's form will be determined 
later. We divide the original Hamiltonian @ into H = H° + H 1 , where H° contains the first 
two terms of H and H 1 the last term. Then the unitary transformation can be proceed order 
by order, 

H' = H° + H 1 + [S, H°] + [S, H 1 } + i[5, [S, H ]] + 0(q 3 ). (7) 
The first-order terms in H' are 



H' + iS,^] = *J29(k + q,k)(bt q + b q )ci +qtS c k>s 

W ^ + Q ' + 9 ' ^^-Q + b Q) C t+q,s C K 



q,k,s 

+-L E 9(k + q,k) 5(k + q, k)(bl q - b q )(e k - e k+q )ct +qjS c kjS . (8) 



q,k,s 



Note that the ground state \go ) of H°, the non-interacting system, is a direct product of a 
filled fermi-sea |FS ) and a phonon vacuum state |ph, ): 



|<7o) = |FS)|ph,0 
Applying the first-order terms on \go ) we get 



(H 1 + [S,H°])\g 



1 



t J 



1 k+q,s 



Ch, 



q,k,s 



l-5(k + q,k)[l 



Cfc — e k+q 



UJ„ 



(9) 



bo), (io) 



,,t 



since 6 g |ph, ) = 0. As the band is half-filled the Fermi energy ep = 0. Thus c k+q s Cfc jS |FS ) / 
only if e k+q > and e k < 0. So, we have 



(H 1 + [S,H°])\g )=0, 



(11) 



if we choose 



5(k + q,k) = 1/(1 + \e k+q - e k \/u q ). (12) 

This is nothing but making the matrix element of H 1 + [S, H°] between \go ) and the lowest- 
lying excited states vanishing. Thus the first-order terms, which are not exactly canceled after 
the transformation, are related to the higher-lying excited states and should be irrelevant 



under renormalization \21\. The second-order terms in H' can be collected as follows: 
[S,H'] + i[S,[S,H ]] 

_ _■_ E E & + *y + i, «) s(k+qM2 _ S{k , + „, k , }] 



q,k,s q' ,k' 
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x(&-g - b q )(bl ql + b q >){c[ +qs c kljS 8 Kkl+ql - c\, +q ,^c k)S 5 kl)k+q ) 
1 \r- \r- g(k + q,k)g(k' + q',k') , , , 

zi\ — — ' 



q,k,sq',k' ~q~q 

1 V- g{k + q,k)g(k' -q,k') , , f f 

~JV^ ^ w <$(fc + g,fc)[2-$(A: - 9, A; )\c k+qs c k , s c k ,_ qs ,c k ,$.2>) 

q,k,sk',s' 1 

5 k > jk+q is the Kronecker 5 symbol. All terms of higher order than a 2 will be omitted in the 
following treatment. 

Second, for the dimerized state, owing to the spontaneous lattice distortion, the neigh- 
boring atoms move in opposite directions. To take into account the static phonon-staggered 
ordering, we make a displacement transformation to H' 

H" = exp{R)H' exp(-i?). (14) 

Here the generator of the displacement operator exp(i?) is 



R = -J2(-l) l u 0] J^(b}-b l ). (15) 

If the ground state of H is \g ), then the ground state of H" is \g' ): \g ) = e~ s e~ R \g' ). 
We assume that for \g' ) the electrons and phonons can be decoupled: \g' ) k, \e) |ph, ), 
where |e ) is the ground state for electrons. After averaging H" over the phonon vacuum 
state we get an effective Hamiltonian for the electrons, 

H cS = (ph,0|F"|ph,0) 

= 2NKul + J2 E o ( k )c{s c k,s + E * A o (k)c^ s c k _^ s 

k k 

"^EE ^ + = ^ S(k + q, k) [2 - 5(k' - q, «)]ci^ t8 c k J„^cm 

q,k,sk',s' 1 



where 



1 ^ g(k',k)g(k ,k') xn/ 

k> ^k'-k 

A (k) = 4au smk[l-6(k-TT,k)]. (18) 



E (k) = e k - ± ]P I ' W, k)S(k, k'){e k - ey ), (17) 



We find by means of the variational principle that the dimerized lattice displacement ordering 
parameter is 

a 



Uq 



KN , 

k>0,s 



E ismk[l-5(k-ir,k)]/e (cJ )S c fc _ W)8 - i^ 8 c k ,s) 



(19) 
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The last term in H e g is a four-fermion interaction. As we are dealing with a one-dimensional 
system, how to treat the four-fermion interaction is a difficult problem. 

Note that in the adiabatic limit, where w f = 0, one has 5(k' , k) = 0, and H e ^ goes back 
to the adiabatic mean-field Hamiltonian, 

H cS (uj n = 0) = 2NKu 2 + ^2 e k <y k jC k ,s + Z ^ au o sinA;(4 )4[ c fc _ 7r)S - c[_ w >s c k}S ). (20) 

k,s k>0,s 

It can be diagonalized exactly by means of the Bogoliubov transformation. The energy-gap 
parameter and the phonon-staggered ordering parameter in this limit can be obtained as 

A(k) = 4au smk, (21) 



m P 



±J2(-iy( Ul )= Uo . (22) 
I 

Thus, we have the relation 

A ad = Aam p , (23) 

same as the results of previous works g, |3|, [f{| in the adiabatic limit. 

On the other hand, in the antiadiabatic limit, where 0% — > oo, we have no = 0, 5(k', k) = 1, 
and H e g becomes 

rr i \ \ 1 g(k + q,k)g(k' - g,k') f | 

k,s q,k,sk',s' ' 

Returning to the real space, this Hamiltonian is 

#eff(^7r -> oo) = -iXX c Ifi, s c M + c L c m,s) 

2 

Z ( c L c m,s + c l+i,s c l,s)(cls' c l+hs' + 4+l, s '<V)- (25) 

i,s,s' 

For the spinless case, the on-site electron-electron interaction disappears because of the Pauli 
principle and this leads the Hamiltonian to be 

2 

H eS (u n oo) = -t ^(cjq+i + cJ +1 q) + ^- ^(cJqcJ +1 q+i - cjq), (26) 
I l 

same as the exact effective Hamiltonian obtained by writing the phonon part of the SSH 
model as a functional integral and integrating out the phonon degrees of freedom IS]. In 
this case, the electrons experience an effective repulsive force between nearest neighbors, and 
for a 2 /2Kt = 1 the system undergoes a transition to a charge-density- wave state. This 
also is equivalent to the antiferromagnetic XXZ spin chain model (via a Jordan- Wigner 
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transformation) with the phase transition taking place at the isotropic point [ 33 1 . Thus our 
effective Hamiltonian works well in both the cj^ = and the uj^ — > oo limits. 

Now the total Hamiltonian can be written as H = Hq + Hi, where H\ includes the terms 
which are zero after being averaged over the phonon vacuum state, and 

# = E^( & k+^)+^ff- ( 27 ) 

Through above two unitary transformations, the Hamiltonian is divided into two parts: Hq 
and H\. We believe that the Hq contains the main physics of the coupling system. This 
is supported by: (a) H c s [Eq. fllq)] works well for the both limits ll> w = and u> n — > oo, 
even if the electron-phonon coupling is strong; (b) the effect of the first order terms in H" 
is eliminated in the lowest order of perturbation by introducing a function 5(k',k); (c) the 
results for the general uj„ in the spinless case are compared quite well [^] with those of the 
Monte Carlo simulations IS]. 



IV Analytical Approach 

Note that the four-fermion interaction term in Eq. ( |i~6| ) goes to zero as uj n — > [see Eq. 
(|20|)]. Therefore, in the case of small u> n , this term can be treated as a perturbation, and the 
unperturbed Hamiltonian is 

H% = 2NKul + £ E (k)c{ s c k>s + J2 iA (k)cl )S c k ^ >s . (28) 

k^s k , s 

The four-fermion term can be re-written as 

H , s _ i g (k + W v - «■ «) s{h + m _ 5{k , _ „ n 

iV g,k,sk',s' U <i 
X ( C t+q,s C k,s — Ck+q-n,s C k-TT,s)(cl,_ qtS ,Ck> t s' ~ Cft -q-ir ,s' C *' '-T ',«') 

1 g{k + q,k- ir)g(k' - q,k' - tt) , , 

Z^ Z^ *(* + Q,k-ir)[2- 5{k -q,k - tt)] 

iV q,k,sk',s' 

X ( C i+g,s C fe-7r,s C I'-(/-7r,s' C fc',s' + C jfc+ 9 -7r,s C M C fc'-g,s' C fc'-7r,s' ) 

1 + 9, fc - ^OfO' - q, k' - tt) , , 

~lvF + q,k-ir)[2-5{k -q,k - tt)] 

7V q,k,sk',s' 

X ( C l+q-Tr,s C k,s c t'-q~iT,s> C k',s' + 4+g,s C fe-7r,s C I'-g,s' C fe'-7r,s' ) • (29) 

In these terms we have the constraints 

k> 0, fc' > 0, fe + q > 0, fc' - q > 0. 
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We can distinguish between different physical processes. The first term in Eq. ( |29| ) is the 
forward scattering one, the second the backward one, and the last the Umklapp one. We 
use the Green's function method to implement the perturbation treatment and it is more 
convenient to work within a two-component representation, 

*k,s = ( ° k,S ), k>0. (30) 



Ck- 



7T,S 



Thus we have 



*l, s ^*fc,s = Cks C k~n,s + A-K,s C k,s > ( 31 ) 

*L io »*M = cLck-^s ~ 4 c k , s 



and the Hamiltonian becomes 



H° eS = 2NKul + E (k)^{ s a z ^ k , s - ]T A (k)^{ s a y ^ k , a , (32) 



and 



rr' 1 sr ST 9(k + q,k)g(k' -q,k') . 

H cs = -^2^2^ 6(k + q,k)[2-S(k -q,k)\ 

q,k,sk',s' " 
^k+ q ,s^kM'- q ,s'^k', S ' 

, 1 9{k + g, k - Ti)g{k' - q, k' - tt) 

+^7 2^ 2^ s \ k + q,k-TT)[2-5( y k -q,k - tt)] 

1 ^ ^ 9(k + q,k-ir)g{k' -q,k' -tt) , , 

x (*i +q ^ x *kM>- q ,s^k>,s> + *I+ 4 , a ^*M*k-,,^v**',«')- (33) 

= x > Vi z ) 1S the Pauli matrix. The matrix Green's function is defined as (the temperature 
Green's function is used and at the end let T — > 0) 



G s , s >(k,T) = -(T T # M (r)*J ,(())> 



= -^5I ex P( _iu; "' r ) G ' s ^'( fc ' a; ™)- ( 34 ) 

P n 

The Dyson equation is 

G s ^(k,uj n ) = G° ss ,(k,u n ) + G° s>a (k,uj n )'E* a> j3(k,L0 n )Gp >s ,(k,uj n ), (35) 

where 

G° s ,(k,u n ) = 6 S;S > {iuj n - E (k)a z + A (k)a y }' 1 , (36) 
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is the unperturbed Green's function. The self-energy S* p(k,u n ) can be calculated by the 
perturbation theory, 

KA k ^n) = ^ E E ^ k ) I 2 - 5 ( k > k')]{G%(k>, Um ) + a z G%(k', co m )a z } 

i r ^ g(k,k' -ir)g(k',k-ir) 

^ k>>0 m ^'-fe-Tr 

x{zo- ?/ G° i( g(A;',cj m )zcj ?/ - (TjG^^',^)^} 

E y 9(k,k-n)g(k',k'-7r) m k _^ + ^ fe/ _ ^ ) _ ^ fc _ k , 
k' m ^ 



x ^T r [ia y G° a f3 (k' , uj m )]ia y + T r {a x G° aJ3 (k' ,uj m )]a x } 



In the theoretical analysis we have taken into account the fact that only the Umklapp scatter- 
ing terms affect the gap, and the forward and backward scattering terms contribute nothing 
to the "charge" gap |34|, |3l|. From Eq. ([37]) one can get that S* g(k,u n ) is irrelative to uj n , 
therefore the spectrum structure of G S;S i(k,u n ) should be 

G SjS < (k, io n ) = S SjS i {iuj n - E(k)o z + t±(k)oyY X ■ (38) 

From G SjS /(k,uj n ) the elementary excitation spectrum in the gapped state can be derived 



W(k) = \J E 2 (k) + A 2 (k). (39) 
The renormalized band function and the gap function are 

E(k) = E,{k)-^Y, 2 4-™ 2 {^)^ k 'M2-K k 'M =ML ,(40) 
N y K V 2 7 ^E 2 (k') + A 2 (k>) 

A(k) = Aauo sin/c[c - d5(k - ir, k)). (41) 

Here 

c = 1 + 1 V ^d(k-7T,k)smk ; A ° (fc) (42) 

N kUs K 2au ^ E 2 {k) + A 2 (k) 

d = 1-4 y ^[l-6(k-7C,k)]smk . A ° (fc) (43) 

N kUs K 2an ov / 'El(k) + A 2 (k) 

The equation to determine uq is 
These are basic equations in our theory. 
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The phonon-staggered ordering parameter 

1 



N 



E(-!)W 



/ 



If u% = we have 6(k',k) = and c = 1, Eq. (41) becomes the same as that in the 
adiabatic theory A a d = A(tt/2) = Aauo. Comparing Eq. (fi"l"|) with the energy gap in 
the adiabatic case, A = 4cmo> we have the gap in the nonadiabatic case, A = A(-7r/2) = 
4auo[c — d]. This is the true gap in the excitation spectrum. 

Fig. 1 shows the elementary excitation spectrum W{k) as function of the wave vector 
in the case of a 2 /Kt = 0.336 for different phonon frequencies uj^/t = 0.001, 0.01, 0.1, and 
0.2. The adiabatic limit result is also shown in dash-dot-dot line for comparison. Near the 
Fermi surface, the change of phonon frequency strongly affects the elementary excitation 
spectrum. The dimerization gap decreases as the phonon frequency increases and changes 
very sensitively to the change of phonon frequency, especially when phonon frequency is 
small. The dimerization gap is much more reduced by the quantum lattice fluctuations than 
other order parameters, for instance, the phonon-staggered ordering parameter. The reason 
for this will be discussed later. In the mean field approximation the Peierls distortion opens 
a gap 2Amf and Amf = 4am p . Our results indicate that this relation holds only in the 
adiabatic limit. The effect of the quantum lattice fluctuations on the dimerization gap can 
be seen clearly also in the mean electron occupation number 

i>u,>=i-!|, (46> 

Fig. 2 shows the mean electron occupation number as function of the wave vector in the 
cases of a 2 /Kt = 0.2 with u> n /t = 0, and a 2 /Kt = 0.4 with u> w /t = and 0.2. The more 
the phonon frequency or the electron-phonon coupling increases the more the Fermi surface 
is smeared. The strongest effect of the quantum phonon fluctuation at the Fermi surface 
leads to the strongest modulation of the dimerization gap. When a = 0, the system becomes 
one for free electrons and there is no long-range order, and consequently, as shown in dash- 
dot-dot line, the system has a clear Fermi surface at k = tt/2. When a > 0, the effective 
electron-electron interactions induced by the electron-phonon coupling correlate the electrons 
and, thereby, the excitation spectrum of the system becomes one of quasi-particle excitations 
superseding the single electron excitations in a = case. 
The density-of-states of fermions p(co) is 

P(") = ^E^-y^ + A^) 

k,s 
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k=f(u), 



(47) 



where, = f{oj) is the inverse function of u = \/E 2 (k) + A' 2 (k). Fig. 3 shows the calculated 
DOS of fermions for electron-phonon coupling a 2 /Kt = 0.336 with uj^/t = 0.001,0.01, and 
0.04, respectively. One can see that a nonzero DOS stars from the gap edge and, for small 
values of uj n , there is a peak above the gap edge with a significant tail between it and the 
true gap edge. As 0% increases, the peak height of the DOS decreases, and the peak broadens 
and moves to lower photon energy which implies that the dimerization gap becomes narrower 
and the Fermi surface becomes more smear. The inverse-square-root singularity at the gap 
edge in the adiabatic case |3[ Q| disappears. For comparison, the adiabatic (uj n = 0) result 
is also shown in the dash-dot line which has an inverse-square-root singularity at gap edge. 
A previous analytic treatment on the modification of DOS due to lattice fluctuation showed 
the move of the DOS peak, but there existed still the singularity [ 10 1 . Our results consist 



with the spectrum measurements in both the peak position and the line shape of DOS. 
The mean phonon occupation number (MPON) is 



1 „ /lt , . 2K 2 1 ^ 4au . A(fc) 1 ^ ^ 4a 2 5 2 (k, k') 2 k + k' 

(48) 

Fig. 4 shows the dependences of MPON on the phonon frequency for a 2 /Kt = 0.336, 0.50 and 
0.81, respectively. MPON decreases as the electron-phonon coupling decreases or the phonon 
frequency increases. For small lo t , MPON changes very sensitively to the phonon frequency, 
and when uj n = 0, MPON is divergent. Therefore, in point of the MPON, adiabatic approach 
is corresponding to the divergence of MPON, while for realistic materials, MPON is finite or 
even small. 



V Optical Absorption Coefficient 

The optical absorption has remained an important source of information about the excitation 
spectrum of (CH)^. |3^], and the measurements of it's optical properties have provided much 
of the experimental basis for theoretical studies |p6[| . In the pristine sample the absorption 
curve peaks at about 2 eV and the peak is followed by a broad tail extending into the gap. To 
appreciate the nonadiabatic effects on the optical excitation spectrum of (CH)^, we further 
calculate the optical absorption coefficient under nonadiabatic circumstance by adopting the 
SSH model. The optical absorption coefficient a(to) can be expressed by the retarded Green's 



function [37| as follows: 

a(uj) = — — ImK R (u), (49) 
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where K is defined as 

K R {u) = 

Here j is the current operator 

j 



- iijt dt( g \\j(o)m-mm}\9)- 



iev F Y,{ci,s c W,s ~ c l+i,s c iy 

2,8 

J sinkW ks a z ^ k>s , 

k>0,s 



(50) 



(51) 



vf = 2at,J = 2evF, and j{t) = exp(iHt)j exp(—iHt) is the form of j in the Heisenberg 
representation. Because the averaging of H% over the phonon vacuum state is zero, in the 
ground state at zero temperature H\ can be neglected. Then, by using the approximately 
decoupling \g') |#q), the ground state of Hq, and H « Hq in the calculation 



li(O)iWb) = ( 5 '|[e(^) je -^)]e^[e^) je -^)]e-^| 5 ') 



(52) 



we can get 



■p 



fc>0,s 



w - 2W(fc) + i0+ w + 2W(k) - i0- 



sin 2 fc - |p £ 9{k '' k } 9 ^ k,) 5 2 (k\ k) sin A-(sin /■• - sin // : 



k-k' 



A 2 (k) 
W 2 {k) 



J 2 ^ 



fc>0,i 



1 



-J, s Vw- Wfe-ifc' - W{k) - W(k') + i0+ oj + oj k _ k , + W(k) + W(fcO - i0+ 



g(k',k)g(k,k') 



UJ 



5 (k ,k) (sin k — sin k 



l\2 



k-k' 



1 — signA: 



,E{k)E{k') + A(fc)A(fc') 



(53) 



Thus, the optical absorption coefficient 



sin 
4J 2 



fc>0,s 
2 



fc 2 V 4a 



k + k' 



-sign sin 



k — k' 



uj n sin A; 



K + 4t|sin^|) 2 



AT2 



8a 2 



k>0,s k' 
a k -j- k 



1 — signAr 



,E(k)E(k') + A(fe)A(fc / ) 



W r (fc)PF(A; / ) 



X ^ COS -2" 



k-k' 



sm • 



K+4i|sin^|) 2 ' 



(54) 



13 



and the to— integrated spectral-weight function 

S(u) = r a(w')du/. (55) 



In Fig. 5, we compare our calculation (solid line) with the observed optical absorption 
(solid circles) [|l3| of trans-(CH) x polymer chain. The optical absorption coefficient a{u) is 
normalized to its peak value cx(u p ) and the photon energy to is in unit of its peak value ui p for 
convenience to comparison. Here we use the input parameters a = 4.1 eV/A, K = 21 eV/A 2 , 
t = 2.5 eV, same as those of Su, Schrieffer and Heeger [||, and uj^/t = 0.0652. One can see 
that the agreement between the experiment result and our calculation is quite good. The 
relationships between the changes of the optical absorption shape as well as its peak position 
and the different phonon frequencies can be seen clearly in Fig. 6. The parameter values used 
are: a 2 /Kt = 0.35, K = 21 eV/A 2 , and t = 2.5 eV with u) n /t = 0.002 (solid line) and 0.04 
(dash line). The calculated true dimerization gaps relating to these two sets of parameter 
values are marked by the arrows. The solid circles denote the experimental results obtained 
from analysis of the inelastic electron-scattering data One can see that the spectrum 
broadens but the peak height decreases and shifts to lower photon energy as w T increases. 
The inverse-square-root singularity at the gap edge in the adiabatic case || 38 1 and in the 
other methods |K| disappears because of the nonadiabatic effect. For finite uj n , there exists 
a significant tail below the peak. This figure shows again the excellent consistency of the 
optical absorption coefficient between our calculation and the experiment result in not only 
the peak position but also the line shape. It should be pointed out that in our calculation 
of the optical absorption coefficient, the choice of the parameter value set for the solid line 
requires a = 4.28 eV/Abeening a little bit larger than a = 4.1 eV/Aused in Ref. [jij, which 
implies that in our theory, due to the nonadiabatic effect, the larger electron-phonon coupling 
is required than that in adiabatic approach to obtain the same optical excitation energy. In 
addition, in experiments the measurement of optical absorption spectrum could be affected 
by various factors, such as the impurity in samples, the finite measurement temperature, 
and the finite measurement resolution. These may make the measured absorption spectrum 
to be broadened and lead to the slight discrepancy between the experiment result and our 
calculation. 

Taking the excellent agreement of our results with the experiment observations on the 
optical absorption coefficient as the check and verification of the effectiveness of our analyt- 
ical method, we further study the difference between the dimerization gap and the optical 
excitation response energy of the trans- (CH) X polymer chain within the SSH model. The 
experiments of many quasi-one-dimensional materials showed that the energy gap deduced 
from the absorption edge is smaller than the activation energy of the dc conductivity fjSfl, 



which, as was pointed out [14, 15], can not be explained in methods assuming a static lattice 
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distortion. A Monte Carlo calculation also suggested that the optical gap can be much larger 
than the dimerization gap Our calculated peak position of optical absorption coefficient 
u; p and the true dimerization gap 2 A = 2A(7r/2) as functions of phonon frequency to n are 
illustrated in Fig. 7. One can see that for the finite phonon frequency, the peak position 
of optical absorption coefficient is not directly corresponding to the true dimerization gap. 
The dimerization energy gap is smaller than the optical gap. It is the most notable that 
there is a discontinuous drop in the dimerization gap once the phonon frequency changes, 
no matter how small it is, from zero to finite, though at the adiabatic limit the dimerization 
gap 2A(-7r/2) = top. After the drop the dimerization gap and the peak position of the optical 
absorption coefficient decrease as the phonon frequency increases. The dimerization gap is 
much more reduced by the quantum lattice fluctuations than the optical gap is. The dimer- 
ization gap is the value of Eq. (41) at the Fermi point k = tt/2, where the quantum lattice 
fluctuations have the strongest effect, while the optical absorption coefficient is the integral 
[see Eq. (|54])] over the all Brillouin zone and the effect of the fluctuations is gentled. We 
note that in our theory, in mathematical viewpoint, the difference of the dimerization gap 
between the uj n = and co n > cases mainly comes from the functional form of the gap [see 
Eq. flU])]. Comparing it with that in the adiabatic limit, one can see that the subgap states 
come from the quantum lattice fluctuations, i.e., the second term in the square bracket of 
Eq. (|4l|). Substituting u> w = into Eq. (|||), we obtain the optical absorption coefficient in 
the adiabatic limit 

gj2 

a(u>) = — -g (4crao) 2 sin 4 k 





-l 


dk_ 





(56) 

k=cf>(u)) 



where, k = <f>(uS) is the inverse function of u = 2ye| + (4auo sin k) 2 . Obviously, k = tt/2 is 
a singular point of Eq. ([56|) , and, at this point, u> p = 2(4aito) = 2A. Thus we get the same 
relation as in the adiabatic approach methods assuming a static lattice distortion, but our 
result indicates that this relation holds only in the adiabatic limit. In nonadiabatic case, in 
our theory, the logarithmic singularity in the integration of Eq. (JH3) in the adiabatic case is 
removed by the factor 1 — 5{k — tt, k) as long as w T is finite. The effects of quantum lattice 
fluctuations on the dimerization gap and on the optical gap are essentially different, especially 
when io n is small. Though some works ascribed the difference between the dimerization 
gap and the optical excitation response energy to the electron-electron interactions [13], our 
calculation reveals that the effects of the quantum phonon can lead the dimerization gap 
to be much smaller than the optical gap in electron-phonon coupling systems. Contrary 
to the statement mentioned in the beginning that the quantum lattice fluctuations cause a 
slightly small reduction of the dimerization order parameter [18, 2C], our analysis shows 
that the dimerization gap being much more reduced behaves very different from other order 
parameters. 
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VI Conclusions 



We have studied the dimerization gap and the optical absorption coefficient of the Su- 
Schrieffer-Heeger model by means of developing an analytical nonadiabatic approach. By 
investigating quantitatively the effects of quantum phonon fluctuations on the gap order and 
the optical responses in this system, we have shown that the dimerization gap is much more 
reduced by the quantum lattice fluctuations than the optical absorption coefficient is. The 
calculated optical absorption coefficient and the density of states do not have the inverse- 
square-root singularity, but have a peak above the gap edge and there exists a significant tail 
below the peak. The peak of optical absorption spectrum is not directly corresponding to 
the dimerized gap. Our results of the optical absorption coefficient agree well with those of 
the experiments in both the shape and the peak position of the optical absorption spectrum. 
When the phonon frequency uj n > 0, there still exists a static dimerization of the lattice but 
the quantum lattice fluctuations (the zero-point lattice motion) play a very important role. 
Taking the SSH model as a typical example for electron-phonon coupling systems in this 
study, we show that the strongest effects of the quantum lattice fluctuations at the Fermi 
surface lead the dimerization gap to be much smaller than the optical gap in electron-phonon 
coupling systems. The difference between the true dimerization gap and the peak position of 
optical absorption spectrum is obtained quantitatively. 
This work was supported by NSF of China. 
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Figure Caption 



Fig. 1 The fermionic excitation spectrum W(k) as function of the wave vector in the 
case of a 2 /Kt = 0.336 for different phonon frequencies uj w /t = 0.001, 0.01, 0.1, and 0.2. The 
adiabatic limit result is also shown in dash-dot-dot line for comparison. 

Fig. 2 The mean electron occupation number as function of the wave vector in the cases 
of a 2 /Kt = 0.2 with u n /t = 0, and a 2 /Kt = 0.4 with u w /t = and 0.2. The result for a = 
is shown in dash-dot-dot line. 

Fig. 3 The calculated DOS of fermions for electron-phonon coupling a 2 /Kt = 0.336 with 
uj w /t = 0.001,0.01, and 0.04, respectively. For comparison, the adiabatic (u, = 0) result is 
also shown in the dash-dot line. 

Fig. 4 The dependences of MPON on the phonon frequency for a 2 /Kt = 0.336, 0.50 and 
0.81, respectively. 

Fig. 5 The comparison of our calculation (solid line) with the observed optical absorp- 
tion (solid circles) of trans-(CU.) x polymer chain. The optical absorption coefficient a(u) is 
normalized to its peak value a(uj p ) and the photon energy oj is in unit of its peak value u p ) 
for convenience to comparison. The input parameters are a = 4.1 eV/A, K = 21 eV/A 2 , 
t = 2.5 eV, and u w /t = 0.0652. 

Fig. 6 The relationships between the changes of the optical absorption shape as well 
as its peak position and the different phonon frequencies. The parameter values used are: 
a 2 /Kt = 0.35, K = 21 eV/A 2 , and t = 2.5 eV with u n /t = 0.002 (solid line) and 0.04 (dash 
line). The calculated true dimerization gaps relating to these two sets of parameter values 
are marked by the arrows. The solid circles denote the experimental results obtained from 
analysis of the inelastic electron-scattering data. 

Fig. 7 The calculated peak position of optical absorption coefficient oo v and the true 
dimerization gap 2 A = 2A(-7r/2) as functions of phonon frequency u n . 
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